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Abstract. 

The magnetic phases of hexagonal perovskites RMn03 (R=Ho, Er, Tm, Yb, Sc, 
Y) are analysed using group theory and the Landau theory of phase transitions. The 
competition between various magnetic order parameters is discussed in the context 
of antiferromagnetic interactions. A phenomenological model based on four one- 
dimensional magnetic order parameters is developed and studied numerically. It is 
shown that coupling of the various order parameters leads to a complex magnetic 
field-temperature phase diagram and the results are compared to experiment. 

PACS numbers: 75.10._b 75.10.Hk 
1. Introduction 

Hexagonal perovskites RM11O3 (R= Ho, Er, Tm, Yb, Lu, Y or Sc) belong to an unusual 
class of materials known as "multiferroics" , which display simultaneously electric and 
magnetic ordering. Most of the hexagonal manganites are ferroelectric below a very 
high temperature (T c « 900 K), and order magnetically at a lower temperature (T N « 
100 K). A complex phase diagram involving different magnetic order parameters has 
been investigated using second harmonic generation [21 C3 Q], neutron scattering 0j 
and heat capacity measurements Perhaps the most intriguing and technologically 
promising development is the recent observation of the magnetoelectric effect in the low 
temperature region of the phase diagram of HoMn03 jHllZj- There are also indications 
of strong magneto-elastic coupling in HoMnOa jHj. 

The manganese and rare earth spins in RMn03 are nearly geometrically frustrated 
because of a slightly imperfect triangular lattice structure, however, complications often 
associated with frustration, such as spin liquid behaviour, are absent here. Each 
magnetic phase is associated with a well-defined non-collinear spin structure within the 
hexagonal plane or antiferromagnetic ordering along the c-axis. The various magnetic 
phases have very similar antiferromagnetic (AF) interaction energies, due to the almost 
perfect triangular lattice, which results in a close competition between them. 
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In each RMn0 3 , the magnetic phase '5 2 ' appears in zero magnetic field at the 
temperature T N) while the phase L develops in magnetic fields of the order of few 
Tesla. These two phases are associated with ordering of Mn 3+ spins within the hexagonal 
plane and are separated by a broad region of hysteresis. Additional phases appear at 
low temperatures in HoMn03 [9j. A sharp Mn spin-reorientation transition occurs at 
Tsr ~ 33 K, which results in the appearance of a third magnetic phase Q B\ [10J. Below 
33 K, a new intermediate phase, which exhibits the magneto-electric effect, has been 
found in the region between the B\ and B 2 phases jHlE]- Moreover, a fourth magnetic 
phase l A\ has been observed below 4 K due to ordering of Ho 3+ spins jT2j. 

In Section II we describe the spin structures associated with each magnetic phase, 
and the antiferromagnetic competition between them. Section III we present a Landau 
model ^31 which describes the four magnetic phases seen experimentally and numerical 
simulations of phase diagrams derived from the model. In Section IV we discuss our 
results. 

2. Crystal and Magnetic Structure 

Below T c , hexagonal RMn03 has the space group symmetry P63cm (#185, C\ v ) [HJ. 
The crystal structure of RMnOa is shown in Figure QJl). With six copies of the chemical 
formula per unit cell, the Mn 3+ ions occupy the (6c) positions, and form triangular 
lattices on the z and z + 1/2 planes. The (6c) positions are (x,Q,z) and equivalent, 
where x ^ a/3 and z — 0. x = a/3 yields a perfect triangular lattice. The rare earth 
ions occupy (2a) and (4b) positions, which are (0, 0, z) (and equivalent, with z = 0.22c) 
and (1/3, 2/3, z) (and equivalent, with z = 0.27c) respectively. A perfect triangular 
lattice is formed when the two z-parameters are equal. The triangular lattices formed 
by the Mn and rare earth ions are shown in Figure El 

The point group Cq v has four one-dimensional irreducible representations (IR's), 
Ai, A 2 , Bi and B 2 , and two two-dimensional IR's Ei and E 2 . The same notation is used 
to label magnetic representations, where characters of £ -l' in the character table of Cq v 
indicate a combination of the point group element with time reversal ^Hj- For the ID 
representations, the same names are given to the corresponding magnetic phases. In 
the presence of magnetic order parameters, the magnetic space groups are P63cm (Ai), 
P6 3 cm (A 2 ), P6 3 cm (B x ), and P6 3 cm (B 2 ). 

The spin configurations of rare-earth and Mn ions may be classified according to the 
magnetic representations by which they transform. All configurations which transform 
according to ID representations are listed in Table [IJ where the spin subscripts refer 
to the atom numbers shown in Figure QJd). The remaining degrees of freedom are 
accommodated by the 2D representations of Cq v ] but so far, there is no evidence 
that phases corresponding to 2D OP's appear in the phase diagram, so the 2D spin 
configurations are not included here. 

We now discuss general features of the H-T phase diagrams of RMnOs by 
considering the relative strength of the nearest-neighbour antiferromagnetic interaction 
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J 1 rtl*rJ 




(a) 



(b) 



Figure 1. a) Atomic positions in a single hexagonal primitive cell of RMnOs. b) 
Numbered Mn ions at the (6c) positions (top) and rare earth ions at the (4b) and (2a) 
positions (bottom). 



for each spin configuration in Table [IJ The antiferromagnetic interaction is defined as 

/=ij£&-4 (i) 

Z <ij> 

where the sum is over nearest neighbours and Si is a spin operator. The parameter J 
depends on the distance between nearest neighbours. 

The first four configurations for Mn ions listed in Table and illustrated in Figure 
confine the spins to the hexagonal plane. The B 2 phase that appears below T c , 
the A 2 phase found in magnetic fields, and the B\ phase found at lower temperature 
in HoMnOa are due to these configurations. In order to compare the AF interaction 
strengths for the different configurations, first, sets of nearest neighbours should be 
separated into three cases, according to symmetry. The first case is the set of co-planar 
nearest neighbours ({1,2,3} and {4,5,6}, numbered as in Figure QJ))). In this case, 
the AF interaction is the same for all four configurations, / = —Js 2 . The second and 
third cases involve non-coplanar pairs. The second case pairs ions on opposite sides of 
the hexagonal primitive cell ({1,4}, {2,5}, {3,6}) and the third case is the remaining 
pairs. The second case favours A\ and A 2 equally, while the third case favours B\ and 
B 2 equally, with / = ± Js 2 in both cases. The distance between partners in each pair for 
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(a) (b) 

Figure 2. Triangular lattices formed by a) Mn ions (viewed parallel to the c-axis) and 
b) rare earth ions (viewed at an angle with respect to the c-axis). 



R (2a) 


A 2 


Slz + S2 Z 




Bi 


Slz — S2 Z 


R (4b) 


Ai 


Slz + S2 Z — Ssz — S^z 




A 2 


Slz + S2 Z + Ssz + S^z 




Bi 


Slz ~ S2 Z — Ssz + S^z 




B 2 


Slz ~ S2 Z + Ssz — S/^z 


Mn (6c) 


Ai 


— (Six ~ S/ix) + \ (S2x + Ss x — S$ x — Sq x ) 






+ ^Y( — S2y + Ssy + Sty — SQy) 




A 2 


— (Siy — Sty) + *y(S2 X — Ss x — S$ x + Sq x ) 






~\~\(S2y + Sty — Sty — Sty) 




Bi 


— (Siy + Sty) + ^y(S2 X — Ssx + S$ x — Sq x ) 






+ 2(^2y + Ssy + Sty + Sty) 




B 2 


— (Six + S± x ) + \ (S2x + Ss x + Ssx + Sq x ) 






+ ^Y( — S2y + Sty — Sty + Sty) 




A 2 


Slz + S2Z + Ssz + S^z + S$ z + Sfo 




Bi 


Slz + $2^ + — S^z — S$ z — Sq z 



Table 1. ID spin configurations classified by irreducible representation. The rare earth 
spin configurations are shown in Figure 0] and the first four Mn (6c) configurations are 
shown in Figure 01 

the second and third cases is exactly the same if the Mn position parameter x is exactly 
1/3. Deviations away from x — 1/3 will favour either the A phases or the B phases. 
However, the observed behaviour is the opposite of what could be expected. At higher 
temperatures (75 K) x = 0.338(1) [4j, which brings the second case pairs closer together 
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and favours the A phases, but the B 2 phase is observed. At the lowest temperature (1.5 
K) x = 0.330(1) the B phases are favoured but the A\ phase is observed. The subtle 
competition between all four phases which results in the dominance of the B 2 phase at 
high temperatures is most likely to be resolved by the inclusion of other interactions, 
such as next-nearest-neighbour, or interactions with the rare-earth ions. 



In HoMn03, additional phases B\ followed by A\ appear as the temperature is 
further lowered. The B\ phase is associated both with in-plane Mn moments and Ho 
ordering along the c-axis. The (2a) and (4b) holmium ions are almost co-planar. The 
B\ phase is AF along the c-axis for (2a) and (4b) ions, but inside the hexagonal plane, 
the (4b) ions are aligned ferromagnetically with each other, and antiferromagnetically 
with the (2a) ions. This is shown in Figure Ek). 

At the lowest temperatures, the A\ configuration emerges, which has AF ordering 
of the (4b) spins in in the planes and along the c-axis, as shown in Figured). Then, 
there can be no AF arrangement with respect to the (2a) positions (because there is no 
A\ configuration for them). The phase transition from B\ to A\ is first-order, as is the 
case between all transitions involving different order parameters, therefore hysteresis is 
anticipated, both on phenomenological grounds, and also because of the persistence of 
the B\ configuration on the (2a) positions. 

The total in-plane AF interaction energy for all Ho ions is approximately the same 
for all three configurations shown in Figure 01 In all cases / = — 2Js 2 , but J differs 
according to the distance between ions. If the (2a) and (4b) position parameters z are 
equal then the three phases are degenerate, otherwise A\ and B 2 have lower energy. 





Figure 3. Manganese spin configurations in the hexagonal plane. These correspond 
to the first four (6c) configurations listed in Tabled (after [J]). 
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(a) (b) (c) 

Figure 4. Rare earth spin configurations for a) Ai, b) B\ and c) B 2 configurations 
for (2a) and (4b) positions listed in Tabled 

The co-linear AF interactions favour the B\ configuration, with / = — 3Js 2 (versus 
I = -2Js 2 for A x and I = 2Js 2 for B 2 ). 

The A 2 phase will always dominate at high enough magnetic fields, since 
it transforms in the same way as the applied field, and couples linearly in the 
phenomenological sense. Microscopically, Zeeman coupling to Mn or Ho spins induces 
ferromagnetic order associated with the A 2 phase. 

3. Landau Model and Phase Diagram 

The order parameters of the phases Ai, A 2 , Bi and B 2 are denoted by 771, ?y 2 , Vs an d 
774, respectively. The minimal Landau model which describes the A 2 and B 2 phases, 
observed in all RMn03, is 

F = OL 2 r}\ + /? 2 ry 2 + a A rjl + far\\ + j^vhl ~ H z {piT] 2 + p 2 rf 2 + p^rjl) (2) 

where where a i) fa and 7^ and pi are phenomenological coupling constants and H z is 
the magnetic field parallel to the c-axis. The coefficients are temperature dependent, 
ai = cii(T — T i ) ) where Ti is the temperature limit of stability for each phase (which for 
convenience, we call the "transition temperature"), and fa > is required for stability. 
<^4 changes sign at T N « 100K in all RMn03. 

In zero applied field, the model allows for four different phases: (0, 0) (the parent 
phase), which may be connected to either (772, 0) (A 2 -phase), (0,774) (B 2 -phase) or 
(^72? ^4) (mixed phase) by second order phase transitions. These are found by solving 
the set of coupled equations dF/drji = 0, subject to the minimisation conditions 
{d 2 FjdT)l) > and (d 2 F/dr) 2 )(d 2 F/dr) 2 ) - (d 2 F / dr) 2 dr] 4 ) 2 > 0. The model also allows 
for the coexistence of two or more different phases by hysteresis. 
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The mixed phase (772,774) can be a minimum of F only when 4/? 2 /?4 > 7I4. Its 
existence is noi the result of hysteresis. Anomalies in the c-axis magnetisation at the B 2 
phase boundary [Sj are evidence that B 2 and A 2 are coupled (i.e. 724 ^ 0). In general, 
r]2 grows linearly with applied field but it may still be subject to a transition in the 
sense that a change in sign of a 2 will increase the number of minima of the Landau 
functional. 

Additional order parameters 771 and 773 are required to describe the additional phases 
observed in HoMn03. Additional terms in the free energy include those obtained by 
replacing in 7?f by 77^ and 773, as well as terms of the form 771772773774 and H^iT]^. 

The Landau model describing Ai, A 2 , Bi and B 2 phases is 

F(Vu V2, ^4) = ^(772, 774) + ai7?i + fcrit + a 3 r]l + /? 3 7?3 + 7l277i7?2 

+ 7l377i773 + 7l477i774 + 723^3 + 734^4 + 7*71*72*73*74 

- H z (p 4 r) 2 r)l + p 5 77277i + PqViVsVi) (3) 

This model may be solved exactly in zero applied field, and it is found that the allowed 
phases are (0, 0, 0, 0) (the parent phase), (771, 0, 0, 0) etc. (Aj or Bj) and (771, 772, 0, 0) etc. 
(mixed phases involving two order parameters). In addition, two or more phases may 
co-exist due to hysteresis. 

The models J2J) and © were analysed by varying the temperature and field and 
searching for the minima of F numerically. Typically, several minima were present, 
the deepest corresponding to the true ground state and the rest to metastable states 
observed as hysteresis. Each set of parameters a ) /?, 7 and p yields a different phase 
diagram; these parameters were varied to find the best match to phase diagrams obtained 
in experiments [IJ El E21 QB| • 

Figure El shows numerical simulations of the phase diagrams for RMn03, modeled 
by ©• In all four diagrams, the onset of the B 2 phase for H = is determined by setting 
T 4 = 80 K. Below this temperature two minima corresponding to non-zero 774 occur in 
(J2J). The A 2 phase, which dominates the right side (high field region) of all diagrams, 
corresponds to a Landau functional which has only one minimum that is shifted away 
from 772 = because of the term ^772 . Thus the free energy resembles the parent 
phase but the symmetry is A 2 . In Figures Eh), c) and d) hysteresis occurs, indicated 
by lighter coloured areas near the black phase boundary line. Here the free energy has 
three minima, one corresponding to the field-shifted parent phase, and two to the B 2 
phase. In Figures Et), b) and d) we have ^ 24 < 4:f3 2 f3 4:) so a pair of minima for 772 7^ 
never occurs. However, in FigureEb), 7I4 > 4/?2/?4, so two shallow minima for 772 co-exist 
with two for 774 in the bottom left corner of the phase diagram. Figure Et>) shows the 
correct arrangement of phases on the phase diagram compared to experiments, but fails 
to simulate correctly the curvature of the phase boundary. Fig ure EtO , which includes 
the non-linear (in OP) field-dependent terms in 0, is a significant improvement. 

Figure El shows numerical simulations of the phase diagram of HoMn03, modeled 
by ©• The upper part of the diagrams, showing the B 2 and A 2 phases, is similar to 
RMn03, in that the B 2 phase appears below the transition temperature T 4 = 75 K, 
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(c) (d) 

Figure 5. Numerical simulations of the phase diagram of RMnOs. The black line 
is the boundary between the A2 and B2 phases. Hysteresis occurs in the lighter- 
coloured regions. In all four diagrams T4 = 80 K, T2 = 10 K, 02 = 10, ^4 = 100 
and pi = 12, while 724 = 60, 150, 300 and 150 in (a), (b), (c) and (d) respectively. 
The parameters ai (o^ = a$(T — Ti)) were scaled to unity. Nonlinear field-dependent 
coefficients p2 = ps = 5 are introduced in diagram (d). 

and the A 2 phase is induced by the magnetic field. Non-zero transition temperatures 
Ti were assigned to all other phases (Ti = 52 K, T 2 = 39 K, T 3 = 62 K) below which 
the phases are metastable and hysteresis occurs. True phase transitions occur below the 
transition temperatures, along the lines where the free energies, evaluated for different 
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phases, are equal. Different shades within a single phase represent the presence of other, 
metastable phases. Figure EK) bears a poor resemblance to experiments. Figure Eb), 
which includes non-linear field-dependent terms, is significantly better, but still fails 
to reproduce qualitatively all of the phase boundary around the B\ phase as seen in 
experiments. 




(a) (b) 

Figure 6. Numerical simulations of the phase diagram of HoMnOs. Black lines 
separate the phases, while different shades within a phase represent the presence of 
other meta-stable states. In both diagrams fa = 0.14, fa = 0.08, fa = 0.2, fa = 0.4, 
while 712 = 1, 7i3 = 1, 7i4 = 1, 723 = 6, 724 = 2.9, 734 = 10, pi = 1 in (a) and 
712 = 90, 713 = 714 = 723 = 724 = 734 = 100, pi = 4 in (b). Diagram (b) also includes 
non-linear field-dependent terms p2 = 1, p3 = 2, p4 = 0.05. 



4. Discussion 

In Section El we argued that even in the absence of true geometrical frustration, there is 
not enough information to predict the magnetic ground state without detailed knowledge 
of the AF interaction strength J. In Landau theory, the microscopic model (pQ) is 
replaced by a phenomenological one (|2|3[) . and the parameter J is incorporated into 
the temperature-dependent c^. Landau theory includes all interactions allowed by 
symmetry, and as such is more general than the AF interaction, which is isotropic. 
However, the proliferation of phenomenological constants also inhibits the predictive 
powers of the Landau model. Nevertheless, Landau modeling is useful because it can 
reveal the minimal elements in a theory that are needed to describe the phase diagrams 
of RMn03. In our analysis, we found that a model based on the usual second and fourth 
order terms, and a linear coupling of the order parameter to the magnetic field, does 
not describe well the observed phase diagrams, especially the curvature of the phase 
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boundaries. The inclusion of non-linear (in OP) field-dependent terms is a significant 
improvement. 

The magneto-electric effect is observed in the region between the B\ and B 2 phases. 
Linear coupling of the magnetic and electric fields of the form a^EiBj can only occur 
when both inversion and time reversal symmetry are absent - these are necessary but 
not sufficient conditions. A mixture of B\ and B 2 order parameters (due to hysteresis) 
does not lower the symmetry enough for magneto-electric coupling. However, domain 
walls, which connect different domains of the same phase, have been implicated in the 
observation of the magneto-electric effect [Hj. Thus OP gradient terms, which couple to 
the magnetic field, may significantly alter the free energy landscape, and could possibly 
replace the non-linear field-dependent terms which we introduced. 

In conclusion, we have studied phase diagrams for RMn03, using group theory 
and Landau theory, by including up to fourth order phenomenological couplings 
between order parameters and non-linear coupling to an applied magnetic field. 
Antiferromagnetic competition between magnetic phases, due to a near perfect 
triangular lattice structure, gives rise to a complex phase diagram. Our simulations 
reproduce the main features seen in experimental results, including the general 
arrangement of phases on the diagram, hysteresis effects and the curvature of phase 
boundaries. 
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